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Abstract 

Until April 2007 the MAGIC telescope used a 300 MSamples/s FADC system to 
sample the shaped PMT signals produced by the captured Cherenkov photons of air 
showers. Different algorithms to reconstruct the signal from the read-out samples 
(extractors) have been implemented and are described and compared. Criteria based 
on the obtained charge and time resolution/bias are defined and used to judge the 
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different extractors, by applying them to calibration, cosmic and pedestal signals. 
The achievable charge and time resolution have been derived as functions of the 
incident number of photo-electrons. 

Key words: fast digitization, FADC, digital filter, Cherenkov imaging telescopes, 
7-ray astronomy. 



1 Introduction 



The Major Atmospheric Gamma ray Imaging Cherenkov (MAGIC) telescope 
[1] uses the IACT technique [2] to study the very high energy (VHE, E > 
50 GeV) 7-ray emission from astrophysical sources, at the lowest possible 
energy threshold. The technique uses Cherenkov radiation: A VHE 7-ray en- 
tering the earth's atmosphere initiates a shower (cascade) of electrons and 
positrons, with a particle density maximum about 10 km above sea level (for 
an energy of 1 TeV). The particles in the cascade produce Cherenkov light 
in a cone of about 1° half- angle, which illuminates an area of around 120 m 
radius on the ground. If the MAGIC telescope is located in this area, part of 
the Cherenkov light will be collected by the telescope mirrors and a shower 
image will be projected onto the photomultiplier tube (PMT) camera. The 
Cherenkov photons arrive within a very short time interval of a few nanosec- 
onds at the telescope camera, whose pixels are fast light sensors such as PMTs, 
so that one can trigger on the coincident light signals. The fluctuations of the 
light of the night sky (LONS) cause background noise. This effect is minimized 
by using low exposure times (signal integration times), typically of the order 
of ten nanoseconds. 

To reach the highest sensitivity and the lowest energy threshold, the recorded 
signals have to be accurately reconstructed. Two quantities are of interest: The 
total signal charge and the signal arrival time. The signal charge (the total 
number of photo-electrons released from the photocathode of the PMT) is 
proportional to the total area below the pulse. The sum of the signal charges 
of all camera pixels is a measure of the shower energy. The signal arrival 
time is given by the time difference between the first recorded FADC sample 
and a characteristic position on the pulse shape, like the maximum, the half- 
maximum on the rising edge or the center of gravity of the pulse. 

The timing information may be used to discriminate between pixels whose 
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signals belong to the shower, and pixels which are affected by randomly timed 
background noise. The pixels with a low signal-to-noise ratio are rejected for 
the subsequent image parameterization [4,5]. 

The main background to 7-rays originates from the much more frequent show- 
ers induced by isotropic hadronic cosmic rays. Monte Carlo (MC) based simu- 
lations predict different time structures for 7-ray and hadron induced shower 
images as well as for images of single muons [3, 18,28,31]. This has two con- 
sequences: On the one hand the arrival time structures across the observed 
Cherenkov shower image, from pixel to pixel, depend on the type of the pri- 
mary particle. On the other hand, also the recorded Cherenkov pulse shape 
inside an individual pixel depends on the primary particle. To exploit the 
pulse shape differences, an ultra-fast digitization of the Cherenkov pulses is 
necessary, as is provided by the most recent upgrade of the data acquisition of 
the MAGIC telescope to a 2 GSamples/s FADC system [11,12]. This paper, 
however, deals with the signal reconstruction of the data taken with the initial 
300 MSamples/s FADC system. Because of its limited sampling speed, we do 
not try to exploit the differences in pulse shape here. 

This paper is structured as follows: In section 2 the read-out system of the 
MAGIC telescope is described, and in section 3 the average pulse shapes of cal- 
ibration and cosmic pulses are reconstructed, from data taken with the FADC 
system. These pulse shapes are compared to those implemented in the MC 
simulation program. In section 4 criteria for an optimal signal reconstruction 
are developed. In section 5 the signal reconstruction algorithms and their im- 
plementation in the MAGIC software framework (MARS [7]) are described. 
The performance of the signal extraction algorithms under study is assessed 
by applying them to pedestal, calibration and MC events (sections 6 to 8). 
Section 9 gives the CPU time requirements for the different signal reconstruc- 
tion algorithms. Finally in sections 10 and 11 the results are summarized and 
an outlook is given. 



2 Signal read-out 

Figure 1 shows a sketch of the MAGIC read-out system, including the PMT 
camera, the analog-optical link, the majority trigger logic and flash analog- 
to-digital converters (FADCs). The response of the PMTs to sub-ns input 
light pulses shows a full width at half maximum (FWHM) of 1.0 - 1.2 ns and 
rise and fall times of 600 and 700 ps correspondingly [8] . A transmitter using a 
vertical-cavity surface-emitting laser (VCSEL) diode modulated in amplitude, 
converts the electrical pulse supplied by the PMT into an optical signal. This 
signal is then transferred via optical fibers (162 m long, 50/125 //m diameter) 
to the counting house [9]. After transforming the light back to an electrical 



4 



signal, the original PMT pulse has a FWHM of about 2.2 ns for a single 
photo-electron pulse, and rise and fall times of about 1 ns. 



receiver board 



low gain 




Fig. 1. MAGIC read-out scheme: the analog PMT signals are transferred via an 
analog optical link to the counting house where - after the trigger decision - the 
signals are digitized by a 300 MHz FADCs system and written to the hard disk of a 
data acquisition PC. 

In order to sample this pulse shape with the 300 MSamples/s FADC system, 
the original pulse is electronically shaped by effectively folding it with a func- 
tion of 6 ns FWHM. In order to increase the dynamic range of the read-out, 
the signals are split into two branches, with gains differing by a factor 10. The 
low-gain branch is delayed by 55 ns and both branches are multiplexed and 
read out by one FADC. The switch from high- to low-gain occurs only if the 
high-gain signal exceeds a pre-set threshold, and 55 ns after this happens. Dur- 
ing the subsequent 50 ns the low-gain signal is connected to the output while 
the high-gain signal is blocked. Figure 2 shows the average reconstructed pulse 
shape (generated by a fast pulser, see section 3) as measured by one FADC. 
A more detailed overview about the MAGIC read-out and DAQ system can 
be found in [10]. 
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Fig. 2. Average reconstructed pulse shape from a fast pulse generator showing the 
high gain and the low gain pulse. The FWHM of the high gain pulse is about 6.3 ns 
while the FWHM of the low gain pulse is about 10 ns. 

The following intrinsic characteristics of the MAGIC read-out system are the 
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most important to affect the signal reconstruction: 

Inner and Outer pixels: The MAGIC camera is constructed with two types 
of pixels, inner and outer pixels, with the following differences: 

(1) Size: The outer pixels have an area larger than the inner pixels by a factor 
four [14]. Their area multiplied by photon detection efficiency, however, 
is higher only by a factor 2.6. 

(2) Gain: The camera is flat-fielded in order to yield a similar reconstructed 
charge signal in all pixels, for the same photon illumination intensity. In 
order to achieve this, the gain of the inner pixels has been adjusted to 
about a factor 2.6 higher than the outer ones [17]. This results in a lower 
charge RMS contribution from the light of the night sky (LONS) for the 
outer pixels. 

(3) Delay: Due to the lower high voltage (HV) settings of the outer pixels, 
their signals are delayed by about 1.5 ns with respect to the inner ones. 

Asynchronous trigger: The FADC clock is not synchronized with the trig- 
ger. Therefore the time t re \ between the trigger decision and the first read- 
out sample is uniformly distributed along the range t re \ e [0, TfadcL where 
T FADC = 3.33 ns is the digitization period of the MAGIC 300 MHz FADCs. 
All FADCs run at the same frequency and phase. The 300 MHz clock signal 
is produced at a central place, multiplicated and distributed by equally long 
cables to the individual FADC modules. 

AC coupling: The PMT signals are AC-coupled at various places in the 
signal transmission chain. Thus the contribution of the PMT pulses due to the 
LONS is on average zero. Only the signal RMS depends on the intensity of 
the LONS. In moonless nights, observing an extra-galactic source, an average 
background rate of about 0.13 photo-electrons per nano-second per inner pixel 
has been measured [11]. 

Shaping: As already mentioned above, the optical receiver board shapes the 
pulse with a shaping time of 6 ns FWHM, i.e. much larger than the typical 
intrinsic pulse width. Since the shaping time is larger than the width of a single 
FADC slice, a strong correlation of the noise between neighboring FADC slices 
is expected. 



3 Pulse Shape 

The fact that the signal pulses are sampled asynchronously by the FADCs al- 
lows one to determine the average pulse shape with high accuracy. To do that, 
the signal samples from different recorded pulses are shifted to a common re- 
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constructed arrival time and normalized to a common area/charge. Therefore, 
the precision of the determination of each point along the pulse shape depends 
on the accuracy of the arrival time and charge reconstruction. Possible biases 
in the charge and arrival time reconstruction may introduce systematic errors, 
whose size are unknown at first hand. Figure 2 shows the average signal from 
a fast pulser as reconstructed by the MAGIC read-out system. The relative 
statistical error of the amplitude value of every reconstructed point is well 
below 10 -2 . The pulser generates unipolar pulses of about 2.5 ns FWHM and 
with a preset amplitude. These electrical pulses are fed to the VCSEL trans- 
mitters, and are transmitted using the same analog-optical link as the PMT 
pulses, and are fed into the MAGIC receiver board. 

Figure 3 (left) shows the normalized average pulse shape for the pulse gen- 
erator in the high and in the low gain, respectively. The intrinsic FWHM of 
the generated pulses is 2.5 ns, whereas it is on average 6.3 ns and 10 ns for 
the pulses reconstructed from the high and low gain chains, respectively. The 
broadening of the low gain pulses with respect to that of the high gain is 
due to the limited bandwidth of the passive 55 ns on-board delay line of the 
MAGIC receiver boards. 
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Fig. 3. Left: Average reconstructed high gain and low gain pulse shapes from a 
pulse generator run. The black line corresponds to the pulse shape implemented 
into the MC simulations [16]. Right: Average reconstructed high gain pulse shape 
for calibration runs with green and UV light (see section 8). All pulse shapes are 
normalized to a common arrival time and area. 

Figure 3 (right) shows the normalized average reconstructed pulse shapes for 
green and UV calibration LED pulses [15] (see section 8) as well as that of 
cosmic events. The shapes of the UV calibration and cosmic pulses are quite 
similar. Both have a FWHM of about 6.3 ns. Since air showers from hadronic 
cosmic rays trigger the telescope much more frequently than 7-ray showers, 
the reconstructed pulse shape of the cosmic events corresponds mainly to 
hadron induced showers. The pulse shape from electromagnetic air showers 
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might be slightly different as indicated by MC simulations [3,18]. The pulse 
shape for green calibration LED pulses is wider and has a pronounced tail. 
The difference between the shapes of the calibration LED pulses is not due to 
the LED light color but to different electronics used for the fast LED drivers. 

The reconstructed pulse shapes for generator pulses, cosmic and calibration 
events permit to implement a representative pulse shape in the MC simula- 
tions, see e.g. the full black line in figure 3, left panel. The shape difference 
between the calibration pulses and the cosmic pulses has to be corrected for 
in the calibration procedure [17]. 



4 Criteria for Optimal Signal Extraction 

The goal of the optimal signal reconstruction algorithm is to compute an un- 
biased estimate of the charge and arrival time of the Cherenkov pulse with the 
highest possible resolution. Let us consider a large number of identical signals, 
corresponding to a fixed number of photo-electrons Aph e - By applying a signal 
extraction algorithm, a distribution of estimated signals iVp hc is obtained, see 
also [33] and references therein. Criteria for an optimal signal reconstruction 
algorithm are developed according to [19]. The deviation between true and 
reconstructed value is given by: 

X = N^ e - N phe . (1) 

The distribution of X has a mean B (the Bias of the estimator) and a variance 
V^.The parameter B is also called the Bias of the estimator and RMSE is the 
Root Mean-Squared Error which combines resolution and bias: 

RMSE = V< X 2 > = VV + B 2 . (2) 

Generally, both B and RMSE depend on N p h e and the background fluctua- 
tions BG. In the case of the MAGIC telescope, the background fluctuations 
are due to the electronics noise and the PMT response to the LONS. The 
signals from the latter have the same shape as those from Cherenkov pulses. 
Therefore, those algorithms which search for the highest sum of a number 
of consecutive FADC slices inside a global time window (so-called sliding 
window algorithms) will have a bias. In case of no Cherenkov signal they 
will typically reconstruct the largest noise pulse. Nevertheless, such a sliding 
window algorithm usually has a much smaller variance and in many cases a 
smaller RMSE than the fixed window extractors, which just sum up a 
fixed number of FADC slices. 

The reconstructed charge should be proportional to the total number of photo- 
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electrons in the PMT. This linearity is very important for the reconstruction 
of the shower energy and hence for the measurement of energy spectra from 
astronomical sources. Deviations from linearity may be caused in different 
ways: At very low signals, the signal will be biased towards too high values 
(positive X). At very high signals, the FADC system goes into saturation, and 
the reconstructed signal becomes too low (negative X). Also, any error in the 
inter-calibration between the high and low gain acquisition channels yields an 
effective deviation from linearity. 

Another important feature of an extractor is its robustness, i.e. its stability 
in reconstructing the charge and arrival time for different types of pulses with 
different intrinsic shapes and background levels: 

• Cherenkov signals from 7-rays, hadrons and muons 

• calibration pulses from different LED color pulsers (with different pulse 
shapes, see figure 3 right panel) 

• pulse generator pulses. 

Finally, the extractor has to accurately reconstruct both the high and low gain 
channels. Due to the analog delay line, the low gain pulse is wider. The total 
recorded time window is relatively small, such that parts of the low gain pulse 
may lie outside of the recorded FADC window. 



5 Signal Reconstruction Algorithms 

We have chosen four algorithms for the study of the reconstruction of the sig- 
nal charge and arrival time: Fixed Window, Sliding Window with Amplitude- 
weighted Time, Cubic Spline with Integral or Amplitude Extraction, and 
Digital Filter. For the signal reconstruction algorithms adopted by other air 
Cherenkov telescopes, see e.g. [29-32]. 

5.1 Fixed Window 

This signal extraction algorithm simply adds the pedestal-subtracted FADC 
slice contents of a fixed range (window) of consecutive FADC slices. The win- 
dow has to be chosen large enough to always cover the complete pulse, other- 
wise physical differences in the pulse position with respect to the FADC slice 
numbering would lead to integration of different parts of the pulse. For this 
reason, the fixed window algorithm adds up more noise than the other consid- 
ered signal reconstruction algorithms. Due to the AC-coupling of the read-out 
chain, the reconstructed signals have no bias. 
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In the current implementation, the fixed window algorithm does not calculate 
arrival times. 



5.2 Sliding Window with Amplitude- weighted Time 



This signal extraction algorithm searches for the maximum integral content 
among all possible FADC windows of fixed size contained in a defined time 
range (global window). The arrival time is calculated from the window with 
the highest integral as: 

, _ l^i=i °l L l /ON 

where i denotes the FADC slice index, starting from slice io and running over 
a window of size ws. The Sj are the pedestal-subtracted FADC slice contents 
and the ti are the corresponding times relative to the first recorded FADC 
slice. 



5. 3 Cubic Spline with Integral or Amplitude Extraction 



This signal extraction algorithm interpolates all the pedestal-subtracted FADC 
slice contents of the full read-out window using a cubic spline algorithm, 
adapted from [20]. In a second step, it searches for the position of the maxi- 
mum of the interpolation function. Thereafter, two different estimators of the 
pulse charge are available: 

(1) Amplitude: the value of the spline maximum is taken as reconstructed 
signal. 

(2) Integral: The interpolation function is integrated in a window of fixed 
size, with integration limits fixed with respect to the position of the spline 
maximum. 

The pulse arrival times can also be estimated in two ways: 

(1) Pulse maximum: The position of the spline maximum determines the 
arrival time. 

(2) Pulse Half Maximum: The position of the half maximum at the rising 
edge of the pulse determines the arrival time. 
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5.4 Digital Filter 



The goal of the digital filtering method [21,22] is to optimally reconstruct the 
charge and arrival time of a signal whose shape is known. Thereby, the noise 
contributions to the amplitude and arrival time reconstruction are minimized, 
see also [24]. For the digital filtering method to work properly, two conditions 
have to be satisfied: 

• The normalized signal shape has to be constant. 

• The noise properties must be constant, i.e. the noise is stationary and in- 
dependent of the signal amplitude. 

As the pulse shape in MAGIC is mainly determined by the artificial shaping on 
the optical receiver board, the first assumption holds to a good approximation 
for all pulses with intrinsic signal widths much smaller than the shaping con- 
stant. Also the second assumption is satisfied to a good approximation: Signal 
and noise are independent and the measured pulse is a linear superposition of 
the signal and noise contributions. LONS conditions. 

Let g(t) be the normalized signal shape (e.g. from figure 3), E the signal in- 
tegral (charge) and r the shift between the timing of the physical and the 
considered/probed signals. Then the time dependence of the signal is given 
by y(t) — E ■ g{t — r) + b(t), where b(t) is the time-dependent noise contri- 
bution. For small time shifts r the time dependence can be linearized. Dis- 
crete measurements yi of the signal at times ti (i = l,...,n) have the form 
Hi = E ■ #j — Et ■ (ji + 0{t 2 ) + bi, where g(t) is the time derivative of the signal 
shape, gi = g(ti) and bi = b(ti). 

The correlation of the noise contributions at times ti and tj can be expressed 
by the noise autocorrelation matrix 

B '■ Bij — (hbj) — (bi)(bj) , (4) 

whose elements can be obtained from pedestal data (see section 7). The noise 
autocorrelation matrix is dominated by LONS pulses shaped by 6 ns FWHM. 
The absolute scale of the matrix elements depends on the LONS level. The 
normalized matrix elements may change by about 10% due to varying LONS 
levels in typical observation conditions. The noise auto-correlation in the low 
gain channel cannot be determined from data. The low gain channel read- 
out is only activated in case the high gain signal is above a certain threshold 
resulting in a measurable low gain signal. It has to be retrieved from Monte- 
Carlo studies instead. 

For a given pulse, E and Et can be estimated from the n FADC measurements 
y = (yi,...,y n ) by minimizing the deviation between the measured and the 
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known pulse shape, and taking into account the known noise auto-correlation, 
i.e. minimizing the following expression (in matrix form): 



X 2 (E, Et) = (y-Eg- Erg) 1 B \y -Eg- Erg) + 0(r 2 ) . (5) 

This leads to the following solution: 

t / x ^, 2x , x (g T B- 1 g)B 1 g-(g T B 1 g)B 1 g 

E = w &mp (t rcl )y + 0(r) , w amp (t Tel ) = {gT B -i g){g T B -i g) _ {g T B -i g? ' 

(6) 

t , x / 2n {g T B~ 1 g)B 1 g - (g T B~ 1 g)B 1 g 

Er = w timc (t rcl )y + 0(r) , ™ timc (t rcl ) = {g T B -i g){g r B -i g) _ {g r B -i g) 2 

(7) 

where £ re i is the time difference between the trigger decision and the first read- 
out sample, see section 3. Thus E and Et are given by a weighted sum of the 
discrete measurements yi with the weights for the amplitude, w a m P (^rei), and 
time shift, ty t ime(^rei), plus 0(r 2 ). To reduce 0(r 2 ), the fit can be iterated using 
g(ti = t — r) and the weights tf a m P /timc(Aci + t) [21]. Figure 4 shows examples 
of digital filter weights. 

The expected contributions of the noise to the error of the estimated amplitude 
and timing only depend on the the shape g(t), and the noise auto-correlation 
B. The corresponding analytic expressions can be found in [21]. 



6 Monte Carlo Studies 



Pulses of a specific number of photo-electrons can be simulated by using the 
Monte-Carlo technique to simulate signal pulses and noise (for the MAGIC 
MC simulations, see reference [16]) Moreover, using MC, the same pulse can 
be studied with and without added noise. In the subsequent studies, the Monte 
Carlo simulation was used to determine, for each of the tested extractors, the 
following quantities: The bias and the charge resolution as functions of the 
input signal charge. 

For the subsequent studies, the following settings have been used: 

• The LONS level in the MC simulations has been set to the value determined 
from extra-galactic source observation conditions (0.13 photo-electrons per 
ns, see [11]). 

• The electronics noise has been simulated without any correlations between 
the FADC samples as a Gaussian distribution with a sigma of 1.6 FADC 
counts (corresponding to about 0.2 photo-electrons) per FADC slice, roughly 
at the level measured in data. Note, that in the data the electronic noise 
introduces a correlation between the FADC samples. 
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Fig. 4. Examples of digital filter weights. Top: cosmic pulses, center: UV calibration 
pulses and bottom: blue and green calibration pulses. On the left side, the high 
gain pulse is shown, one the right side, the low gain. Full lines show the normalized 
signal shapes g(t) (multiplied by 5 for better visibility), dashed lines the amplitude 
weights w amp (t), and dotted-dashed lines the time weights wti me (t). F° r the high- 
-gain extraction 4 FADC slices are used and for the low-gain extraction 6 FADC 
slices. 
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• The intrinsic arrival time spread of the photons was set to be 1 ns (FWHM 
of a Gaussian), as expected for 7-ray showers. 

• The conversion factor from photo-electrons to integrated charge over the 
whole pulse was set to 7.8 FADC counts per photo-electron. 

• The relative timing between the trigger and the signal pulse was uniformly 
distributed over 1 FADC slice. 

• The total dynamic range of the entire signal transmission chain was set to 
infinite, thus the detector has been simulated to be completely linear. 



6. 1 Bias 



The signals were simulated with noise and extracted using the different ex- 
tractor algorithms. For all sliding window algorithms the extraction window 
was allowed to move 5 FADC slices, independently of its size. For each signal 
extraction algorithm the average conversion factor between the reconstructed 
charge in FADC counts and the input number of photo-electrons was deter- 
mined separately. The signal reconstruction bias was calculated as a function 
of the simulated number of photo-electrons iV s j m : 

B =< N TCC - N sim > (8) 

Figure 5 shows the results for some tested extractors, with different initial- 
izations. As expected, the fixed window extractor does not show any bias up 
to statistical precision. All other extractors however, do show a bias. Usually, 
the bias vanishes for signals above 5 photo-electrons, except for the sliding 
window. In this latter case, the bias only vanishes for signals above 12 photo- 
electrons. 
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Fig. 5. Charge reconstruction bias as a function of the number of generated pho- 
to-electrons from MC simulations including electronic noise plus LONS. Above 12 
photo-electrons, the bias vanishes for all signal extractors. 
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6.2 Root Mean Square Error 



In order to obtain the precision of a given extractor, we calculated the relative 
RMSE: 

Rel. RMSE = —— ^Vai[N I0C ] + B 2 . (9) 

1 ' sim 

Figure 6 shows the relative RMSE for the high gain and low gain parts sepa- 
rately. Also the square root of the relative variance of the number of simulated 
photo-electrons (\Jl/ N sim ) is shown, which corresponds to the intrinsic fluc- 
tuations of the signal from air showers, following Poissonian statistics. Note, 
that the PMT introduces an additional excess noise [6], which is on average 
18% of the Poissonian fluctuations for the MAGIC PMTs. For all extractors 
the variance of the reconstructed signal is dominated by noise and only slowly 
increases with rising signals due to mis-reconstruction of the signal pulse itself. 
Therefore, the relative RMSE is proportional to 1/iVphe- For small numbers of 
photo-electrons, extractors with small extraction windows or the digital filter 
yield the smallest values of RMSE, but the difference is only important below 
about 5 photo-electrons. Above that value, the curves for all extractors have 
crossed the black line, i.e. they are more precise than the intrinsic fluctuations 
of the signal. This is also true for the entire low gain extraction range. The 
best results are obtained with the digital filter or a spline integrating 1 FADC 
slice. 
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Fig. 6. Relative RMSE as a function of the number of generated photo-electrons 
from MC simulations including fully simulated electronic noise plus LONS. Left: 
high gain, right: low gain. The black dashed line shows the square root of the 
relative variance of the incoming numbers of photo-electrons, note, that the PMT 
introduces an additional excess noise [6]. The best results are obtained with the 
digital filter or a spline integrating 1 FADC slice. 
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7 Pedestal Extraction 



The pedestal is the average FADC count for the signal baseline (no input sig- 
nal). To determine the pedestal setting off-line, dedicated pedestal runs are 
used, during which the MAGIC read-out is triggered randomly. The fluctua- 
tions of the signal baseline are due to electronics noise and LONS fluctuations. 
Thus the pedestal RMS is a measure for the total noise level. 
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Table 1 



The statistical parameters square root of reconstructed signal variance, bias, RMSE and A"pj^, cs - for the tested signal extractors, applied 
to pedestal events. All units are in reconstructed numbers of photo-electrons, statistical uncertainty: about 0.1 photo-electrons. The 
extractors yielding the smallest values for each column are marked in red. 



By applying the signal extractor to pedestal events, the bias B and the RMSE 
for the case of no signal (iVph e = 0) can be determined. Table 1 shows the 
bias, the square root of the variance Var[iV rec ] and the root-mean-square er- 
ror for randomly triggered pedestal events with closed camera and sample 
observations outside the Galactic plane ("extra-galactic LONS", 0.13 photo- 
electrons/ns) as well as within the Galactic plane ("galactic LONS"). In ad- 
dition, table 1 shows the corresponding values from MC simulations. The 
Var[iV rec ] of the MC simulations is slightly lower then in the sample observa- 
tions outside the Galactic plane, although the simulated LONS level in the 
MC simulations has been adjusted to the "extra-galactic LONS". This is in 
part due to neglecting the correlation of the FADC slices from the electronic 
noise, see section 6. Every extractor window had the freedom to move 5 FADC 
slices, i.e. the global window size was fixed to five plus the extractor window 
size. 

One can see that the bias typically decreases and the variance increases with 
increasing sliding window size, except for the digital filter. The extractor with 
the smallest RMSE is the digital filter fitting 4 FADC slices with an RMSE 
of 1.4 and 1.7 photo-electrons for an extra-galactic and a galactic star- field, 
respectively. 

In the so-called image cleaning [5] only the camera pixels above a certain 
charge threshold are used for the image parameterization [4]. The charge 
threshold is adjusted such that the probability of being a noise fluctuation 
does not exceed a certain value. For the sake of comparison, a typical value of 
3 u (0.3% probability) was chosen here and that number approximated with 
the formula: 

^thrcs. ^ B + 3 . ( 1Q) 

iVp|^ cs ' is shown in the 11 th and 15 th column of table 1. Most of the sliding 
window algorithms yield a smaller signal threshold than the fixed window, 
although the former ones have a bias. 

The lowest threshold of only 4.2 photo-electrons for the extra-galactic star- 
field and 5.0 photo-electrons for the galactic star-field is obtained with the 
digital filter fitting 4 FADC slices. This is almost a factor 2 lower than the 
fixed window results. 



8 Calibration 

In this section, tests are described which were performed using light pulses of 
different color, shape and intensity produced by the MAGIC LED calibration 
pulser system [17]. Such a system is able to provide fast light pulses of 2-4 ns 
FWHM, with intensities ranging from 3 to more than 600 photo-electrons in 
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one inner PMT of the MAGIC camera. These pulses can be produced in three 
colors: green, blue and UV. Table 2 lists the available colors and intensities. 



Color 


Wavelength 


Spectral Width 


Min. No. 


Max. No. 


Secondary- 


FWHM 




[nm] 


[nm] 


Phe's 


Phe's 


Pulses 


Pulse [ns] 


Green 


520 


40 


6 


120 


yes 


3-4 


Blue 


460 


30 


6 


600 


yes 


3-4 


UV 


375 


12 


3 


50 


no 


2-3 



Table 2 



The pulser colors available from the calibration system 

Whereas the pulse shape of the UV LEDs is very stable from event to event, 
the green and blue LED pulses can show smaller secondary pulses about 10- 
40 ns after the main pulse. Note, that the UV-pulses are only available in 
intensities which do not saturate the high gain read-out channel. However, 
the brightest combination of (blue) light pulses easily saturates all high gain 
channels of the camera, but does not saturate the low gain read-out. 



8.1 Number of Photo-electrons 

The mean number of photo-electrons < iV p h e > was calculated for a sequence of 
calibration pulses of same intensity, following the excess noise factor method [6] 
and using different signal extractor algorithms. If the signals are extracted 
correctly, < > should be independent of the signal extractor. 

In our case, an additional complication arises from secondary pulses of the 
green and blue colored light pulses, which may introduce a dependence of 
< ^Vphe > on the extraction time- window size (recall figure 3). For the standard 
MAGIC calibration procedure [17] only UV calibration pulses are used. 

Figure 7 shows < iVp hc > for the standard UV calibration pulse. The results 
differ by less than 5%, which results in an additional systematic error to the 
absolute energy scale of the reconstructed events. Note, that the total system- 
atic error of the absolute energy scale was estimated to be 16% [34] . A small 
increase of < iVphc > for an increasing window size can be observed. This may 
be due to the intrinsic time structure of the calibration pulse. 

The peak-to-peak variation of the conversion factor between FADC counts and 
number of photo-electrons for the different intensities is below 10% [25] for all 
extractors. The corresponding non-linearity is due to the intrinsic non-linearity 
of the MAGIC signal chain and a possible non-linear signal extraction. 



19 



50 r 

45: 



40 r 

35 lm 



A 


30 


a 






25 


V 


20 




15 




10 



Inner Pixels, rel.RMS=0.045 



CM 4 <S CO 

3 3 3 3 



5 5 5 S 

O) O) O) OJ 



5-S 

: x 



3 (0 (0 (0 (0 

0.0)0)0)0) 
g <D O 0) OJ 

< £ £ E £ 

(U <D <D <D O 



extractor 



u-tntntntntntnintntn 



100 
80 
t 60 

a 

40 
20 




Outer Pixels, rel.RMS=0.042 



CM 



CD CO 



• > » J J ? : 51 ! 5 

o "o "o "o "o = ™ ™ ^ 

= = = = 

13 .E .E .E .E c c c c c 
SJ-o-o-o-u:=:=:=:=:= 
.2 = = = = a. a. a. a. a. 
u-cntocncocnaitnujcn 

extractor 



Fig. 7. < iVphe > from the standard calibration pulse, emitted by 10 UV LEDs, 
reconstructed with different signal extractors. Left: inner pixels, right: outer pixels. 
The statistical errors are smaller than the marker size. 

8.2 Robustness Tests 



Possible variations of the pulse shape may degrade the signal extraction 
quality of the MAGIC data. Variations of the pulse form have a physical 
reason: average Cherenkov pulses from hadronic showers are usually broader 
than those from electromagnetic cascades. Additionally there are differences 
between the pulse form of calibration pulses and those of cosmic pulses. These 
variations affect mainly those signal extractors which integrate only parts of a 
pulse or perform fits to a sample pulse form. In order to quantify the magnitude 
of the effect, table 3 lists the fraction of the pulse which is contained in typical 
time windows around the pulse maximum, for various pulse forms. While the 
amplitude extraction or integration of only 1 FADC slice around the maximum 
yield differences as high as 10% (cosmic - UV), the error is reduced to about 
3% if four FADC slices are being integrated. Further deviations, characterized 
by the blue LED calibration pulse or the MC pulse, yield an even stronger 
discrepancy. 

The digital filtering method assumes a constant signal shape to compute the 
weight functions. In fact, all pixels are assumed to have the same average 
signal shape and the same weights are used for all pixels. In order to test 
the robustness of the digital filtering method with respect to deviations of the 
actual pulse shape from the assumed pulse shape, the standard UV calibration 
pulse was extracted using different weight functions (computed for UV and 
blue calibration pulses, Cherenkov pulses and the low gain). The results are 
displayed in figure 8 showing variations of about 8% in the reconstructed 
signal for typical pulse form variations (blue, UV and cosmic weights) and 3% 
in < iVphe > after calibration using the same weights. Note, photons from a 
7-ray shower (from a weighted spectrum) arrive within (2-2.5) ns [3]. 
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Window Size 






High Gain 




Low Gain 


(FADC slices 


MC 


Cosmic Calib. UV 


Calib. Blue 


Cosmic 


Calib. Blue 


around maximum) 






(percentage of complete pulse integral) 




Amplitude 


5.4 


5.0 


4.5 


4.1 


3.4 


2.5 


1 slice 


54 


50 


46 


41 


35 


27 


2 slices 


78 


76 


71 


66 


60 


48 


4 slices 


97 


98 


95 


89 


90 


82 



Table 3 



Pulse form dependency of integration windows: Shown is the fraction of the signal 
(in percent of the complete pulse integral), contained in different time windows 
around the pulse maximum for different pulse shapes. In the case of the first line 
"Amplitude" , the signal amplitude has simply been divided by the complete pulse 
integral (arbitrary units). 

In conclusion, an event-to-event variation of the pulse shape may cause a 
charge reconstruction error of up to 10% for all signal extractors which do 
not integrate the entire pulse. The size of this error decreases with increasing 
integration window size. For the digital filter this event-to-event variation of 
the pulse shape may lead to an error of the reconstructed signal of up to 8%. 
A systematic difference in the pulse shape causes an error of up to 3% (after 
calibration) in the case of the digital filter. 

8.3 Time Resolutions 

The calibration light pulses can be used to test the time resolution of sig- 
nal extractors. Thereby, the arrival time difference St is measured for every 
channel, with respect to a reference channel: 

SU =U — t re f , (11) 

where ti denotes the reconstructed arrival time of pixel number i and i re f the 
reconstructed arrival time of a reference pixel. Using a calibration run of a 
fixed number of calibration pulses, the mean and RMS of the distribution of 
Sti for a given pixel can be computed. The RMS is a measure of the com- 
bined time resolutions of pixel i and the reference pixel. Assuming that the 
photomultipliers and read-out channels are of the same kind, an approximate 
time spread of pixel i is obtained from the sigma of a Gaussian fit to the 
distribution of the time differences Stf 

At, « a(Sti)/V2 . (12) 
Figure 9 shows the obtained average time resolutions < At, > as a function of 



21 



Inner Pixels, rel.RMS (high gain)=0.076 



to 

> 

3 





0) 


in 


to 


*t 


10 






*- 




2: 


.n 


£ 




13) 




'o> 


'(I) 


■5 






3 




a> 


> 




3 


3 


m 


CO 



in 

'5 

o 

E 

in 
o 
O 



in 
£ 

C) 

'5 
3 
o 

I 
in 
o 
O 



0) 
■* 

D) 
'(I) 
3 
C 

'w 

D) 

o 



in 
10 

- 

D) 

'5 

- 

c 

'as 

_ 

3 
o 



Digital Filter extractor 



350 
300 
250 

£200 

v 

150 
100 
50 




Outer Pixels, rel.RMS (high gain)=0.077 



m 

> 
=5 





in 


<n 


m 






10 








£ 




.c 


O) 


D) 


D) 


■<1) 


'<!) 


'O 


3 


3 


3 


0) 


0) 


> 


3 


3 


3 


m 


m 



0) 

■* 
sz 

D) 
'<D 

3 
u 

E 

in 
o 
U 



to 

£ 

'0 

o 
E 

U) 

o 
O 



to 

g 

3 
c 
'5 

O) 

o 



- 

- 
: 
- 

- 
: 
- 

: 
- 



Digital Filter extractor 



Inner Pixels, rel.RMS (high gain)=0.019 



£ 

- 

3 



■ 


■ 


■ 






w 


0) 




ID 


10 




*- 




£ 


f 


£ 




C7) 


O) 


■(1> 


'o> 


■5 


3 


3 


3 


01 


o> 


> 


3 


3 


3 


CO 


CO 



in 
£ 

O) 

'5 
3 
o 

E 

</> 
o 
o 



in 
10 

£ 

D) 
'5 

3 

u 
I 

</) 

o 
O 



in 



in 
10 



.c £ 

3 3 
c c 




Digital Filter extractor 



Digital Filter extractor 



Fig. 8. Mean reconstructed charge in FADC counts (top) and < iV p hc > (bottom) 
from a standard calibration pulse reconstructed with a digital filter using different 
weight functions (computed for UV and blue calibration pulses, Cherenkov pulses 
and the low gain). Left: inner pixels, right: outer pixels. The relative RMS was 
calculated for the first six (high gain) pulse forms. The statistical errors are smaller 
than the marker size. 

< -ZVphc > for various calibration runs taken with different colors and light in- 
tensities for the telescope pointing outside the Galactic plane ( "extra-galactic 
LONS", 0.13 photo-electrons/ns). Three time extractors were used: a Sliding 
Window of 6 FADC slices with amplitude-weighted time, the Cubic Spline 
with the position of the half-maximum at the rising edge of the pulse as ar- 
rival time and the digital filter. Note, that a time resolution of better than 1 ns 
can be obtained for all pulses above a threshold of 5 photo-electrons. For the 
largest signals, a time resolution as good as 200 ps can be obtained. In order 
to understand the exact behavior of the time resolution, we briefly review the 
main contributions: 

(1) The intrinsic arrival time spread of the photons on the PMT: This time 
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Fig. 9. Reconstructed mean arrival time resolution as a function of the ex- 
tracted mean number of photo-electrons for three different extractors: the ampli- 
tude-weighted sliding window with a window size of 6 slices (left), the half-maximum 
searching spline (center), and the digital filter with correct pulse weights over 6 slices 
(right). Error bars denote the spread (RMS) of time resolutions of the investigated 
channels. The marker colors show the applied pulser color, except for the last (green) 
point where all three colors were used. The lines are a fit using equation (17), see 
text for details. The best fit parameters are shown as an inset. 

spread can be estimated roughly by the intrinsic width 5t m of the input 
light pulse. The resulting time resolution is given by: 



At 



6t 



IN 



N, 



(13) 



phc 



The width Stm is about 1 ns for 7-ray pulses, a few ns for hadron pulses, 
for muons a few hundred ps and about 2-4 ns for the calibration pulses. 
(2) The transit time spread (5txTS of the photo-multiplier (the spread of the 
times between the release of an electron from the photo cathode and the 
corresponding signal at the PMT output) which can be of the order of a 
few hundred ps per single photo-electron, depending on the wavelength 
of the incident light. As in the case of the photon arrival time spread, the 
total time spread scales with the inverse of the square root of the number 
of photo-electrons: 



At 



phc 



(3) The reconstruction error due to the background noise and limited extrac- 
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tor resolution: 

5t rcc ■ R/phe . . 

At « rec /P (15) 

JVphe 



where i? = y Var[iVp he ] is the square root of the extractor variance, which 
depends only very weakly on the signal charge. 
(4) A constant offset due to the residual FADC clock jitter between different 
channels or the MC simulation time steps. 

At « St . (16) 

In total, the time spread can be expressed as: 



AT 



n t, 



2 



\ iVp he jV£ he 



+ T$. (17) 



where Ti contains the contributions of Stj^ and St^TS, the parameter T 2 con- 
tains the contribution of <5t rec and To the offset St®. 

The measured time resolutions in figure 9 were fitted by equation (17). The 
low fit probabilities are partly due to the systematic differences in the intrinsic 
pulse shapes of the different color LED light pulses. Nevertheless, all calibra- 
tion colors had to be included in the fit to cover the full intensity range. In 
general, the time resolutions for the UV pulses are systematically better than 
those for the other colors. This can be attributed to the fact that the UV 
pulses have a smaller intrinsic pulse width [25] and the UV LEDs are very 
stable from event to event, whereas the green and blue LED pulses can show 
secondary pulses about 10-40 ns after the main pulse (see section 8.1), which 
influence the reconstructed pulse arrival time. 

There are clear differences between the studied time extractors, especially the 
sliding window extractor yields poorer resolutions. This is in part due to the 
fact that in the chosen sliding window algorithm the bias of the time recon- 
struction with respect to the relative timing (phase) of the pulses with respect 
to the free running FADC clock has not been corrected for. The parameters Ti 
and To should in principle be independent of the time extraction algorithm. 
Nevertheless, T\ is larger for the sliding window algorithm than for the spline 
interpolation and the digital filter. This is in part due to the (anti)-correlation 
between the reconstructed charge and arrival time for the former extractor, 
see equation (3). 

From the measured time resolution for calibration pulses one can estimate 
the expected time resolution for cosmic pulses. The only important difference 
between calibration and cosmic pulses are different arrival time spreads of the 
photons on the PMT camera. The time spread of the photons on the PMT for 
cosmic pulses is smaller than for blue/green calibration pulses, but about the 



21 



same as for the UV calibration pulses (see the widths of the pulses in figure 
3b). Therefore, the timing resolution for cosmic pulses is at least at the level 
of the timing resolution determined from the calibration pulses. The timing 
resolution for cosmic pulses is conservatively estimated to: 



AT 



\ 



4.5 ns 2 20 ns 2 



iVi 



phc 



phe 



0.04 ns 2 . 



(1* 



For signal charges above 10 photo-electrons the time resolution is below 830 ps. 
For signals of 100 photo-electrons the time resolution may be as good as 300 ps. 



9 CPU Requirements 



The speed of different extractor algorithms (the number of reconstructed 
events per unit time) was measured on an Intel Pentium IV, 2.4 GHz CPU 
machine. Table 4 shows the average results whereby the individual measure- 
ments could easily differ by about 20% from one try to another (using the 
same extractor). The numbers in this list have to be compared with the event 
reading and decompression speed (400 events/s). Every signal extractor being 
faster than this reference number does not limit the total event reconstruction 
speed. Only some of the integrating spline extractor configurations lie below 
this limit and would need to be optimized further. 



Name 


Events/s. 




(CPU) 


Fixed Window 8 slices 


3200-4000 


Sliding Window 6 slices 


1000-1300 


Spline Amplitude 


700-1000 


Spline Integral 1 si. 


300-500 


Digital Filter 


700-900 



Table 4 

The extraction speed measured for different signal extractor configurations. Note, 
that the fixed window does not calculate the arrival time. 



Note, that the calibration light pulses illuminate directly the camera, whereas the 
cosmic light pulses are reflected by the MAGIC mirror system. The MAGIC mirrors 
have been built in a parabolic shape, and are thus isochronous. Nevertheless, they 
have been staggered in a chess-board manner [26] with an offset of about 10 cm. 
This introduces an additional contribution of about 700 ps width to the intrinsic 
arrival times spread of the Cherenkov photons. 
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10 Results and Discussion 



The results based on the investigations discussed above are summarized in ta- 
ble 5. Note, that there is no absolute basis for criteria to separate acceptable 
from non-acceptable properties of signal extraction algorithms. In the follow- 
ing the arbitrarily chosen criteria to compare the extractors are motivated: 

• The extractor should yield on average the true number of photo-electrons 
and should not deviate by more than 10% in case of slight modifications of 
the pulse shape. These deviations directly effect the determination of the 
absolute energy scale of the reconstructed events. Note, that the dominant 
systematic error to the absolute energy scale is currently the photon detec- 
tion efficiency (10-12%) [34]. This requirement excludes extractors which 
integrate only a small portion of the pulse, especially the amplitude sensing 
cubic spline extractor. 

• The extractor must yield a stable low gain pulse extraction. This means 
that apart from being robust against modifications of the pulse shape, the 
extractor has to reconstruct on average the true signal charge also in case 
of variations of the pulse position within the recorded FADC samples. This 
criterion excludes the fixed window extractor since arrival time jitters may 
exceed the time window between the tail of the high gain pulse and the 
beginning of the low gain pulse. 

• The RMSE of the reconstructed charge for the case of no signal should not 
exceed 2 photo-electrons (an arbitrarily chosen threshold) for dark night ob- 
servations and the RMSE of the reconstructed charge for air shower signals 
should never exceed the intrinsic Poissonian signal fluctuations plus excess 
noise above 5 photo-electrons. Camera pixels with a signal below 5 photo- 
electrons are usually rejected for the image parameterization [4,5]. This 
low-energy analysis condition discards the large sliding windows and the 
fixed window extractor. It is not critical for high-energy analyses, however. 

• For analyses close to the energy threshold, an extractor should have a small 
or negligible bias, discarding again the amplitude sensing cubic spline ex- 
tractor. 

• The time resolution should not be worse than 2 ns at a signal strength of 
10 photo-electrons. Note, that this condition allowed us to require a time 
coincidence of 3.3 ns between neighboring pixels to reject noise signals in the 
image cleaning and thus allow to measure differential energy spectra down 
to 60 GeV [34]. All fixed window and all simple sliding window extractors 
are excluded by this condition. 

• The needed CPU-time should not exceed the one required for reading the 
data into memory and writing it to disk. Unless further effort is made to 
speed up the integrating spline, it is excluded if used with a large integration 
window. 
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Extractor 


robust- 


robust- 


RMSE 


bias 


time 


Speed 


Configuration 


ness 


ness 






spread 






pulse 


pulse 












form 


form 












high gain 


low gain 










Fixed Window 8 si. 


OK 


NO 


NO 


BEST 


NO 


BEST 


Sliding Window 2 si. 


NO 


NO 


OK 


OK 


NO 


OK 


Sliding Window 4 si. 


OK 


NO 


NO 


OK 


NO 


OK 


Sliding Window 6 si. 


OK 


OK 


NO 


OK 


NO 


OK 


Sliding Window 8 si. 


OK 


BEST 


NO 


OK 


NO 


OK 


Spline Amplitude 


NO 


NO 


OK 


NO 


OK 


OK 


Spline Integral 1 si. 


NO 


NO 


OK 


OK 


BEST 


OK 


Spline Integral 2 si. 


NO 


NO 


OK 


OK 


BEST 


OK 


Spline Integral 4 si. 


OK 


NO 


OK 


OK 


BEST 


NO 


Spline Integral 6 si. 


OK 


OK 


NO 


OK 


BEST 


NO 


Digital Filter 4 si. 


OK 


NO 


BEST 


OK 


OK 


OK 


Digital Filter 6 si. 


OK 


OK 


OK 


OK 


OK 


OK 



Table 5 



The tested characteristics for every extractor. See text for descriptions of the indi- 
vidual columns. OK means, the extractor has passed the test, NO that the extractor 
failed and BEST that the extractor has succeeded a particular test as best of all. 



Table 5 shows which extractors satisfy the above criteria. One can see that 
there is no signal extractor without problems. However, the digital filter fitting 
four FADC slices can always be used for the high gain extraction, and the 
digital filter fitting six FADC slices for the low gain extraction. the mean pulse 
position is not critical. This combination has been chosen as the standard 
extractor for all MAGIC data before April 2007 with the 300 MSamples/s 
FADC system. During a certain period, the pulse position was by mistake 
shifted with respect to the FADC read-out samples. In this case [27] the signal 
was reconstructed by the cubic spline algorithms integrating 1-2 FADC slices. 



If efficiencies at low energies are not critical, i.e. a high analysis energy thresh- 
old without the use of the timing information, the sliding window extractor 
can be used in configurations which cover the entire pulse. This extractor turns 
out to be especially robust and was used for the data analysis in [35]. 
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11 Conclusions and Outlook 

In this paper different algorithms to reconstruct the charge and arrival time 
from the FADC read-out samples of the MAGIC telescope have been devel- 
oped. These algorithms are tested using MC simulations, pedestal and calibra- 
tion events. The achievable charge and arrival time resolutions are determined. 
A digital filter fitting four FADC slices in the high gain channel and six FADC 
slices in the low gain was chosen as the standard signal extraction algorithm. 

Part of the difficulties to find a suitable signal extractor (reflected in table 5) 
stem from the fact that the MAGIC signals are shaped just as long as to cover 
about four FADC slices. This choice, although necessary for a 300 MSamples/s 
FADC read-out, "washes out" the intrinsic pulse form differences between 
7-ray and hadron showers, and thus prevents the analysis from using this 
information in the 7/hadron discrimination, see [36]. On the other side, the 
shaping time is not long enough to safely extract the amplitude from the 
(shaped) signal. 

These problems will be overcome with the full analysis of data taken with 
the new 2 GHz FADC system in MAGIC [11,13]; for first results see [37]. 
This system has been designed to reduce any pulse form deformation to the 
minimum. It can be expected that the individual pulse forms are then directly 
recognized as such, e.g. with a digital filter using two sample pulse forms (a 7- 
ray-like and a hadron-like) and discriminating between the two with the help 
of the calculated \ 2 . These FADCs have a higher dynamic range and do not 
need a separate low gain channel any more. It can be expected that the signal 
extraction will become more robust, besides extracting a wealth of additional 
information about the shower characteristics. 
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